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Abstract 

The size effect in ferroelectrics is treated as a competition between the geometrical symmetry 
of the ferroelectric sample and its crystalline symmetry. The manifestation of this competition 
is shown to be polarization rotation, which is driven by temperature and/or size variations, thus 
providing a qualitative indication of intrinsic finite size effect on ferroelectrics. The concept is 
demonstrated in a simple case of PbTiC>3 nanowires having their axis parallel to [Hl]c direction, 
where the size effect is shown to include polarization rotation and an additional phase transition. 

PACS numbers: 77.80.B-, 77.80.bj 
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Thermodynamic phase transitions in finite systems and relevant size effects have been 
a focal point of interest of experimentalists and theorists for a long time. Nowadays, it is 
not a rarity, that an "ordered" (e.g. ferroelectric) state is documented and manipulated 
in systems having their smallest dimension in the range of a few lattice constants 3|, 
a scale on which the phase transition can really "feel" that the system is finite. At the 
same time, the identification of the nature of a size effect still remains a challenge due to 
the multitude of factors that influence the "ordering" in objects with limited dimension. 
For example, for ferroelectric phase transitions, one can count factors like depolarizing field 
and direct coupling between the polarization and the interface, which can be classified 
as associated with the intrinsic size effects, and factors like size- dependent stoichiometry 
and microstructure, which can be classified as the extrinsic size effects. The experimental 
identification of the nature of size effects in a system is often complicated by the fact that 
most of the theoretical work in the field has been focused on intrinsic scenarios, whereas 
experimentally, extrinsic effects often compete with those intrinsic. On the practical level, 
the problem is further complicated by the fact that modeling of intrinsic size effects shows 
typically a lowering of the transition temperature T c s|-6], a result that is qualitatively 
identical to the expected manifestation of extrinsic size effects. Thus, to prove that the 
phase transition in a given system is prone to an intrinsic size effect, one has to analyze the 
functional size dependence of T c , taking into account also information (typically lacking) on 
the impact of the extrinsic effect on T c . In this context it seems valuable to identify features 
of the intrinsic size effects that would enable their distinction on a qualitative level from 
the extrinsic ones. We show below such a feature of the intrinsic size effect and illustrate it 
using a simple model of ferroelectric nanowires. 

Consider the intrinsic size effect in a sample of a so-called "cubic" ferroelectric, i.e. a 
ferroelectric with a cubic paraelectric phase such as the classical perovskites BaTiC>3 and 
PbTiOa. In the bulk form, the ferroelectric instability in such systems can be classified 
as "isotropic". It is controlled by the matrix of the inverse susceptibilities which is 
proportional to the unitary matrix: 

^ = mop, = a6 - (1) 

where Pi stands for the Cartesian components of the polarization and <3> is the free energy 
density of the ferroelectric. In the Landau theory, a is a temperature-dependent constant 
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changing its sign at the point where the paraelectric phase loses its stability. In the general 
case, the unstable direction with respect to appearance of the spontaneous polarization is 
controlled by the eigen vector of xTj corresponding to the smallest eigen value (which, with 
lowering of the temperature, changes the sign first, from positive to negative). In the case 
of a cubic symmetry, because of the degeneracy of the eigen values of x7/ > this matrix 
cannot govern the directions of the ferroelectric axes of the crystal, which are controlled 
by higher derivatives of the free energy gp^ffkap, and mdP^W^aK • However the 
situation changes essentially when the phase transition of a finite cubic ferroelectric sample 
is addressed. The couplings that control the intrinsic size effect can contribute to the matrix 
Xlj 1 = q^qp. , where the free energy $ comprises now also these couplings. Clearly the 
symmetry of such a contribution is governed by that of the sample or by the superposition 
of the crystalline symmetry and the geometrical symmetry of the sample. This implies that, 
in general, in a finite sample, the ferroelectric instability is not isotropic any more. It will be 
of the easy-plane or easy-axis type, depending on the degree of degeneracy of the smallest 
eigen value of xTj 1 matrix. 

Here of a remarkable interest is the situation where the orientation of the symmetry 
elements of the Xij 1 matrix is different from those of the point group symmetry of the 
crystal. Under such conditions, tracing a phase transition affected by a pronounced size 
effect (as a function of temperature or of a certain sample dimension), one can pass from 
a regime where the correction to the bulk value of Xij 1 * s large, to a regime where this 
correction is negligible (virtually bulk regime). In doing so, one passes from a regime where 
the orientation of the polarization is controlled by the geometrical symmetry of the sample to 
that where it is controlled by the crystalline symmetry. The implication is that, in general, 
the intrinsic size effect on phase transitions in cubic ferroelectrics should be accompanied 
by a rotation of the spontaneous polarization. 

To illustrate this effect, we perform a simple modeling of the ferroelectric phase transition 
in nanowires of PbTiC>3. Let us consider the ferroelectric phase transition in a free-standing 
PbTiC>3 nanowire of a cylindrical shape with its geometrical axis parallel to the [111] di- 
rection (hereafter orientations of crystallographic directions and planes will be given in the 
cubic reference frame of the paraelectric cubic phase). Since, in the bulk from, the ferroelec- 
tric axes of PbTi0 3 are parallel to (100) directions, off- axis components of the polarization 
are expected in these nanowires. We assume that the depolarizing field associated with 
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these components |— d by environments charges. The resmts of experimental flfl 
and theoretical [3, |8| studies suggest that such screening can be very efficient, however it 
is obviously incomplete. The origin of the non- vanishing depolarizing field under screening 
by ambient charges is the seperation of the screening charges from the bound charges (see 
e.g. [3]). The depolarizing energy associated with the incompleteness of such screening is a 
driving force for size effect on the phase transition of this system. Other driving forces for 
this size effect are the short-range coupling of the order parameter with the interface 
and the surface tension^]. In order to provide a transparent illustration of the concept 



formulated above, we retain only the depolarizing effect in our consideration [ll| and incor- 
parate the effect of an incomplete polarization screening in a ferroelectric nanowire within 
the Landau theory framework. 

Addressing the problem of incomplete polarization screening in a short-circuited ferro- 
electric parallel-plate nano-capacitor, the depolarizing effect is described in terms of an 
additional potential jump across the interface, which is proportional to the out-of-plane 
component of the polarization [14j and to the effective screening length A, which character- 
izes the effective thickness of the double electric layer formed by the bound and free charges 

50 = —p . n (2) 

where n is unit vector along the surface normal and Eq is the permittivity of the free space. 
For the perovskite/metal interface, both simple quantum-mechanical 15[ and advanced ab 
initio |4| estimates evaluate A as a fraction of Angstrom. Assuming that for the case of 
screening with environmental charges the values of A are comparable to those given above, 
we conclude that the thickness of the double electric layer formed due to the polarization 
screening in nanowires are much smaller than the typical values of their radii (2-50 nm 
e.g. from js]). Bearing this in mind, we apply Eq.© locally on the side surface of the 
nanowire. Since outside the nanowire the electrostatic potential must be constant, Eq.([2j) 
yields (to within a constant) the boundary conditions for the potential inside the nanowire. 
The solution to the Laplace equation \/ 2 (p = with such boundary conditions leads to a 
uniform depolarizing field, which is antiparallel to the polarization component P± normal 
to the axis of the nanowire: 

^dep = ~^P± (3) 
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where R is the radius of the nanowire. The corresponding depolarizing energy reads 
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The incorporation of this correction into the energy of the system modifies its instability to 
the easy-axis type. In the cubic reference frame, the modified now reads 
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with a = T q^° where C and To are the Curie- Weiss constant and Curie- Weiss temperature, 
respectively. With such input, the system can be described using the following expansion 
for the free energy of mechanically free nanowire 



$ = - 2 \,/ + a n (P x 4 + P 2 4 + P 3 4 ) + a l2 (P X 2 P 2 2 + P 2 2 P 3 2 + Ffl») + a m (Pf + P 2 6 + P 3 6 ) 
+ a 112 [Pi 4 (P| + PI) + P 3 4 (Pi 2 + Pf) + P 2 4 (P! 2 + P 3 2 )] + a^s^PlPl (6) 

with xij 1 coming from Eq.^3]) and the high-order terms taken in the standard form for 
perovskite ferroelectrics. Here, we have used the total polarization P as the order parameter 
in the free energy expansion, though the problems involving depolarizing effect require, 



strictly speaking, the use of the "soft-mode" part of the polarization for this purpose |l6 |. 
However, one can readily show on the lines of Ref . jl6!] that such substitution is legitimate if 
the nanowire radius R is much larger than Ak& where Kb is the background relative dielectric 
permittivity of the ferroelectric. With a typical value of Kb for perovskites of about a few 



units [12||, we expect this requirement to be fulfilled for ferroelectric nanowires. 

For PbTi0 3 wires oriented parallel to [111] direction, the equilibrium states have been 
determined by minimizing the free energy $ with respect to the 3 components of the polar- 



izations vector P. The results of the calculations 18 1 are presented in FigJT] (the absolute 
value of P and its orientation as functions of temperatures T and the wire radius expressed 
in a unit of the effective screening length R/X) and in Figj2] (the phase diagram of the system 
in T - R/X coordinates). It is seen that there are three regimes of behavior of the system, 
illustrating competition between the symmetry of the geometrical sample and its crystalline 
symmetry. 
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FIG. 1: (color online), (a) The amplitude of P s and (b) the orientation of P s indicated by the 
tilting angle from [001] (tetragonal) towards [111] (rhombohedral) in the (110) plane as a function 
of T in a PbTi03 nanowire having its geometrical axis along [111]. The accompanying numbers 
indicate R/X. Note: Curves "2000", "500" and "400" partially overlap with the curve "300" at the 
high-temperature region; The tilting angle of P s for R/X = 300 shows that P s is parallel to [111]. 



(I) For R/X smaller than some 340, it is the sample symmetry which mainly controls 
the system, which behaves as an easy-axis rhombohedral ferroelectric with its polarization 
parallel to the wire axis at all temperatures. 

(II) For R/ X larger than some 7200, it is the crystalline symmetry that mainly controls the 
system. It behaves as multi-axial ferroelectric with 6 possible pseudo-tetragonal orientations 
of the polarization, which is very close to those in bulk PbTiOs. In this regime the size effect 
manifests itself in a weak rotation of the spontaneous polarization in {110} planes passing 
through [111] direction and a shift of the transition temperature Tq from its bulk value down 
to T . 

(III) For intermediate values of R/X, a crossover regime takes place. Here, a first-order 
transition (driven by variations of temperature or/and wire radius) between the rhombo- 
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FIG. 2: The phase diagram of PbTiO"3 nanowires with the geometrical axis along [111]. The 
thick solid line indicates the first-order phase transition and the thin solid line corresponds to the 
second-order phase transition. The dashed lines are the spinodals of the first-order phase transition. 



hedral and pseudo-tetragonal phases takes place. In this case, a jump-like rotation of the 
polarization at the transition and a further substantial rotation in the pseudo-tetragonal 
phase (driven by variations of T and R/X) occur as a manifestation of the competition be- 
tween the geometrical and crystalline symmetries. The spinodals of this transition are also 
shown in Figj2] 

A remarkable feature of the phase diagram shown in Fig{T^, is that the ferroelectric 
phase transition between the paraelectric and the rhombohedral phases is of a second order 
whereas in the bulk material, the transition is of a first order. The origin of this feature can 
be elucidated using the polar plot (Figj3j) of the magnitude of the P -contribution to the free 
energy for a fixed value of P as a function of the orientation of P. It is known that the sign 
of this magnitude controls the order of the transition. It is seen that, for (100) directions, 
the sign is negative, controlling the first-order phase transition in bulk PbTiC>3. At the same 
time, for (111) directions, the sign is positive, explaining the second-order phase transition in 
"thin" nanowires where the depolarizing effect urges the appearance of polarization parallel 
or antiparallel to [111] direction. 

It is interesting to estimate, based on FigJIJ the dimensions of nanowires that are likely 
to manifest a strong size effect. If one assumes that the screening length is 0.1 A, the curves 
representing R/X of 300, 400, 500, and 10000 represent wire diameters of 6, 8, 10, and 
200 nm, respectively. In other words, a nanowire of a realistic 8 nm diameter is predicted 
to show, upon cooling from the paraelectric phase, a second-order phase transition to a 
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FIG. 3: The magnitude of the P 4 -contribution to the free energy for a fixed value of P as a function 
of the orientation of P in the {110} plane in PbTi03. The dashed and solid curves indicate negative 
and positive sign respectively. 



rhombohedral ferroelectric phase, followed by another phase transition (rhombohedral to 
monoclinic) at about 80 °C (FigJTJ). Experimental results are not available at present for 
PbTiC>3 nanowires with growth direction along (111). In the case of Pb (Zrm s Tin.s.O Qa , 
monodomain wires were obtained, having [001] parallel to the nanowire axis 19(. However, 
albeit having a rather large diameter (d^50 nm), polydomain nanowires of KNbC>3 20], 
having (100) direction along the wire axis, showed a low temperature monoclinic phase 
instead of the rhombohedral phase that occurs in bulk KNbC>3, which might be interpreted 
as a fingerprint of an intrinsic size effect in line of the arguments outlined above. 

All in all, the presented model clearly illustrates the qualitative feature of the intrinsic 
size effect on phase transitions in finite samples, i.e. an order parameter rotation occurs in 
the cases where the geometry of the sample lowers the symmetry of the system compared to 
the crystalline symmetry. Despite the simplicity of the model we believe that it provides a 
qualitatively correct picture of the phase transition in PbTiC>3 nanowires of [111] orientation 



211 ] . Keeping in mind a fraction of Angstrom as an estimate for A, the model predicts, for 



realistic dimensions (R lying in the range of several to few tens of nm), essential rotation of 
the spontaneous polarization caused by variations of the temperature and wire radius R. It is 
important to note that short-range polarization/interface coupling and elastic effects related 
to the surface tension could also break the isotropic nature of the ferroelectric instability, 
similarly leading to a rotation of the spontaneous polarization during the phase transition 
in cases where these effects manifest themselves. As the rotation of the spontaneous polar- 
ization entails a rotation of the spontaneous deformation, this qualitative manifestation of 
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the size effect might be detected by XRD analysis. 
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